BABAR-CONF-01/99 
SLAC-PUB-9139 
February, 2002 



(N 

Rare Hadronic and Radiative Penguin B Decays at BABAR 

(N ' 
<D 

Qn ■ 

Stephane Willocq 
Department of Physics 
University of Massachusetts, Amherst, MA 01003, USA 

<N 

, Representing the BABAR Collaboration 

(N 
O 
CI 

o 

">< 

<D . 

JL • Abstract 

i-C ■ We report recent results in the study of rare hadronic and radiative penguin decays of B mesons. These 

results are based on a sample of 23 million BB pairs collected by the BABAR Collaboration at the SLAC 
PEP-II e+e- B Factory. 

rS 

S3 



Contributed to the Proceedings of the 5 KEK Topical Conference: 
Frontiers In Flavor Physics (KEKTC5), 
20-22 Nov 2001, Tsukuba, Ibaraki, Japan. 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



Rare Hadronic and Radiative Penguin B Decays at BABAR 



S. Willocq a representing the BABAR Collaboration 

a Physics Department, University of Massachusetts, Amherst, MA 01003, USA 

We report recent results in the study of rare hadronic and radiative penguin decays of B mesons. These results 
are based on a sample of 23 million BB pairs collected by the BABAR Collaboration at the SLAC PEP-II e + e~ 
B Factory. 



1. MOTIVATION 

A major effort is currently underway to test 
the flavor sector of the Standard Model (SM) in 
general and the origin of CP violation in par- 
ticular. Measurements of the time-dependent 
CP-asymmetries in B° — > J/tpK® decays have 
led to the first observation of CP violation jlj] 
in the B°-B° system and to a measurement of 
the phase (3 of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element V^-the current world av- 
erage is sin 2(3 = 0.79 ± 0.10. The magnitude of 
Vtd is constrained by studies of B° and B® mixing 
and is expected to be determined with an uncer- 
tainty of ~ 5% in the near future. 

The determination of the magnitude and phase 
(7) of the CKM element V u b is expected to be 
much more difficult, due in part to the small rate 
for b — > u transitions (~ 0.5 — 1% compared to the 
rate for b — > c transitions). In parallel with the 
measurement of sin 2(3 with B° — > J/ipK® decays, 
one can in principle determine sin 2a from the 
time-dependent asymmetry in B° — > ir + ir~ de- 
cays (where a = tt— (/?-(- 7)). However, the study 
of charmless two-body B decays has shown that 
gluonic penguin b — > d processes contribute signif- 
icantly to the overall B° — > 7r + 7r~ decay rate and 
have a different weak phase than the dominant 
b — » u tree process, thereby preventing a straight- 
forward interpretation of the measured asymme- 
try. Improved understanding of gluonic penguins 
is thus of very high importance and constitutes 
an important challenge for current experiments. 



Another motivation to study rare hadronic and 
radiative penguin decays is that many of these 
decays proceed via effective flavor-changing neu- 
tral current (FCNC) transitions. Such FCNC 
processes are forbidden at tree level and involve 
heavy particles in loops (e.g. W boson or top 
quark). This implies that the rates are low but 
also that such processes are sensitive to the pres- 
ence of heavy non-SM particles in the loops. 

Besides the measurement of decay rates, one 
can also study CP violation in the decay (a.k.a. 
direct CP violation). In this case, the relevant 
observable is 

r(.B —*•/) + T(B — » /) 

If both penguin and tree decay amplitudes con- 
tribute to the decay, the asymmetry can be ex- 
pressed as 

2|P||r|sinA^sinAg 
cp |p|2 + | T |2 + 2 |p||T|cosA<£ cosAtf ' 1 ' 

where \P\ and |T| represent the magnitudes of 
the penguin and tree decay amplitudes, whereas 
Acf> and AS correspond to the difference between 
the weak and strong phases for the penguin and 
tree amplitudes, respectively. From this expres- 
sion it can be seen that Acp is expected to be 
very small in the SM (less than 1%) for penguin- 
dominated decays like B — > <f>K(*\ K°ir, K*j, 
and thus the observation of large CP asymmetries 
in these modes would clearly signal the presence 
of new physics in the penguin loops. 



2. EXPERIMENTAL DETAILS 

In the following, we report the results of analy- 
ses based on a sample of 23 million BB events 
collected by the BABAR Collaboration at the 
asymmetric-energy PEP-II e + e~ collider. This 
sample corresponds to 20.7 fb _1 of e + e~ annihi- 
lation data recorded at the T(45) resonance. An 
additional (off-resonance) sample corresponding 
to 2.6 fb _1 recorded 40 MeV below the resonance 
is also used to study continuum backgrounds. 
Charged particles are tracked and measured with 
the combination of a five-layer double-sided sili- 
con microstrip detector and a 40-layer drift cham- 
ber immersed in a 1.5 T solenoidal magnetic field. 
Surrounding the tracking systems, a Cherenkov 
ring imaging detector (DIRC) relying on total in- 
ternal reflection of the Cherenkov light produced 
in 4.9 m-long quartz bars provides charged par- 
ticle identification. A Thallium-doped Csl elec- 
tromagnetic calorimeter is used for photon detec- 
tion and electron identification. Finally, an in- 
strumented flux return provides muon and K® 
identification. Further detail about the BABAR 
detector is given elsewhere M . 

Most of the analyses presented here attempt 
to fully reconstruct the B decay. To do so, we 
rely on the well-known beam energies and exploit 
kinematical constraints from energy and momen- 
tum conservation in the e + e~ — > T(4S) — > BB 
process. Two mostly uncorrelated kinematical 
variables are used: the beam-energy substituted 
mass, m E s = \jE^ am - p^ 2 , and the energy dif- 
ference in the e + e~ center-of-mass system (CMS), 
AE = E* B - E£ eam , where £ & * eam and E* B are the 
energies of the beam and the B meson candidate, 
and p|j is the B meson candidate momentum vec- 
tor, all quantities being defined in the CMS. In 
some analyses, a kinematical fit is performed on 
the momenta of the B decay products by impos- 
ing E* B = E'l eam and the corresponding energy- 
constrained mass rriEC is used. The tues resolu- 
tion is dominated by the beam energy spread and 
has a typical value of 2.5 MeV. The AE resolu- 
tion is different for each particular final state and 
ranges between 20 and 70 MeV (the resolution is 
noticeably worse for modes containing photons or 

7T°S). 



With the exception of the B -> K^i+i' anal- 
ysis, the background from BB events is small or 
negligible for all other analyses. The dominant 
background from e + e _ — > qq continuum events 
(q = u, d, s, c) is suppressed using event shape and 
energy flow variables. These variables exploit the 
fact that BB events are much more spherical than 
continuum events due to the small B momentum 
(~ 340 MeV/c) in the CMS. 

Analyses have been optimized and tested us- 
ing signal and background Monte Carlo samples, 
as well as off-resonance data and side-band data. 
The signal region remained hidden during this op- 
timization and validation process. 

3. RADIATIVE PENGUIN DECAYS 

As mentioned earlier, the measurement of B 
and B® oscillation frequencies will soon provide a 
determination of the ratio |Vtd/Vt s | with an un- 
certainty of ~ 5%. Radiative penguin decays are 
also sensitive to the same CKM elements and can 
provide a complementary determination via the 
ratio T(B -► pj)/T{B -> K*~f). 

A first step toward such a measurement is 
the determination of the branching fraction for 
B — > K*j decays. The analysis proceeds with the 
selection of high-energy photons with 1.5 < E 1 < 
4.5 GeV and 2.30 < E* < 2.85 GeV in the labora- 
tory and CMS frames, respectively. Photons con- 
sistent with originating from ir° and rj decays are 
removed. Photons are further required to be iso- 
lated to remove contamination from high-energy 
7T° decays. K* candidates are reconstructed in 
K+n-, K° s tt°, K+tt° and K°tt+ modesQ, making 
use of DIRC information to select charged kaons. 
The dominant backgrounds consist of events with 
photons from initial state radiation or from 7r° 
and r\ decay. These are suppressed by utilizing 
event shape variables, the polar angle of the B 
candidate and the helicity angle of the K* decay. 

B — > K*"f candidates are required to satisfy 
m ES > 5.2 GeV and -0.2 < AE < 0.1 GeV (for 
K+tt- and K a s it+ modes) or -0.225 < AE < 
0.125 GeV (for K° s ir° and K+ir° modes). An 
unbinned likelihood method is used to extract 

1 Charged conjugate states are implied throughout, unless 
noted otherwise. 
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Table 1 

Summary of signal yields and branching fraction 
measurements for radiative penguin decays. The 
first and second uncertainties in the branching 
fractions are statistical and systematic, respec- 
tively. The uncertainty in the yield is statistical 
only. 

Mode Signal yield B (1CT 6 ) 

135.7 ± 13.3 42.4 ±4.1 ±2.2 
14.8 ±5.6 41.0 ± 17.1 ±4.2 

28.1 ±6.6 30.1 ±7.6 ±2.1 

57.6 ± 10.4 55.2 ± 10.7 ±3.8 



-1 



11 \i\ 



< 1.7 (90% C.L. 



K£+£-~^ 0.0 ±0.3 < 0.6 (90% C.L.) 
K*l+t- 0.7 ±1.1 < 2.5 (90% C.L.) 



Figure 1. Distribution of ities for B — > K*^ can- 
didates in four different K* decay channels. 



signal yields from the ities distributions, see 
Fig. [y. The branching fractions for each of the 
four modes are shown in Table |l|. Averaging the 
two K*° and K*+ modes yields B(B° -> K*°j) = 
(42.3±4.0(stat)±2.2(syst)) x 10~ 6 and B{B+ -> 
K*+-y) = (38.3 ± 6.2(stat) ± 2.2(syst)) x 10~ 6 . 
These values, although somewhat lower, agree 
with recent theoretical predictions, given the ~ 
40% uncertainty in the predictions. 

We have also searched for the rare decay B° — > 
77, which is predicted to have a branching frac- 
tion in the range of (0.1 — 2.3) x 10~ 8 within 
the SM. However, this could be enhanced by new 
physics contributing to the loops. Given the cur- 
rent size of the event sample, any observation 
would be clear evidence for physics beyond the 
SM. The analysis proceeds in a way similar to the 
K*^f analysis. One event is found in the signal re- 
gion, defined by the ±2cr ranges \ities — mB° \ < 
0.0078 GeV and \ AE\ < 0.144 GeV, in agreement 
with the expected background of 0.9 event. As 
a result, a limit is set to be B(B° — > 77) < 
1.7x 10~ 6 at the 90% C.L. This represents a factor 
of 20 improvement over the best previous limit. 



We have searched for B + — » K^ + ( + £ with 
K*+ -> K°ir+ and B° -> K^i+i- with K*° -> 
K + tt~. Not only does this analysis need to sup- 
press background from continuum events, it also 
needs to suppress background from BB events. 
Major sources of BB background are events in 
which both B mesons decay semileptonically or 
one decays to J/ipK^ or ip(2S)K^ final states 
with J /if) -> £+£- or V(2S) -> £+£-. It is partic- 
ularly important to suppress B decays to charmo- 
nium since these have characteristics that are sim- 
ilar to the signal. The event yields are extracted 
with an unbinned likelihood fit in the AE-rriES 
plane. The overall efficiencies for the different 
modes range between 6% for K* + fi + fi~ and 18% 
for K + e + e~ . No significant signal is observed but 
the limits (see Table [l]) are close to SM calcula- 
tions. 

4. RARE HADRONIC B DECAYS 

Rare hadronic decays have been searched for in 
two-body, quasi two-body and three-body modes. 
In the case of two-body decays, the following 
topologies are examined: h + h'~, h + ir , K®h + , 
K®ir° and K®K®, where h represents a charged 
pion or kaon. Signal yields are extracted using an 
extended likelihood function taking all relevant 
topologies into account. The likelihood incorpo- 



rates information from event kinematics (AE and 
m^s), event shape and energy flow via a Fisher 
discriminant, as well as mode-specific information 
like the Cherenkov angle for charged hadrons, the 
invariant mass for K® — > ir + ~n~ candidates or the 
helicity angle. This analysis benefits from the ex- 
cellent n/K separation provided by the DIRC. 
The significance of the separation is greater than 
4 a up to laboratory momenta of 3 GeV/c and is 
2.5 (7 at the highest momentum (~ 4.3 GeV/c). 

Signal yields and branching fractions are sum- 
marized in Table ^. The large rate for Kit final 
states as compared to tttt final states confirms the 
previous observation by the CLEO Collaboration. 
Since tree diagrams for B° — > K + ir~ decays are 
Cabibbo-suppressed as compared to B° — > n + n~ 
decays, the larger rate for Ktt final states is a 
clear indication that gluonic penguin diagrams 
contribute significantly to the decay. 

A host of quasi two-body decays were also stud- 
ied. These involve u>, i], if or <p resonances and 
offer another window to study decays mediated by 
either or both tree and penguin decay amplitudes. 
For example, B — ► urn decays are dominated by 
tree decay amplitudes, whereas B — > <j)K de- 
cays are dominated by gluonic penguin decay am- 
plitudes. The later decay mode also holds the 
promise of a sin 2/3 measurement via the time- 
dependent asymmetry in B° — > (f)K®. Such a 
measurement is of interest to test the consistency 
of the CKM picture in two different decay mech- 
anisms, B° — > J/ipK® involving b — ► ccs transi- 
tions and B° — > (f>K® involving b —> sss transi- 
tions. B — ► <j)K and B + — » cf>K* + decays were 
first observed by the BABAR Collaboration. A 
significant signal is also observed in B + — > ujtt + 
decays, see Table § and Fig. § 

Previous studies of B — > rf'^K^ decays re- 
ported a surprizingly large branching fraction. 
With a larger data sample, we confirm higher 
than expected rates for B — > rjK* and B — » r/'K 
decays, see Table | and Fig. |. Since tree con- 
tributions are Cabibbo-suppressed one needs to 
invoke a large enhancement in the penguin con- 
tributions, possibly originating from interference 
between different penguin amplitudes (g* —* uu 
and g* — > ss). 



Table 2 

Summary of signal yields and branching fraction 
measurements for rare hadronic decays. The first 
and second uncertainties in the branching frac- 
tions are statistical and systematic, respectively. 
The uncertainty in the yield is statistical only. 
Mode Signal yield B (1(T 6 ) 
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Figure 2. Distributions of rriEC an d AE for 
(a,b) B+ -> (c,d) B° -> 7/if , and (e,f) 

£> + — > ijjtt + . The shaded area corresponds to 
77' — > ?77r7r. The solid curves represent the result 
of the likelihood function and the dashed curves 
correspond to the background contribution. 



To further explore the large decay rate into 77', 
we have performed a semi-inclusive reconstruc- 
tion of the decay B — ► r/'X s , where X s represents 
a hadronic system with non-zero strangeness. 
The 77' candidates are reconstructed in the t]tt + tt^ 
channel (with 77 — + 77) and are required to have 
momentum p* > 2 GeV/c in the CMS to re- 
move most of the b — ► c — ► rf background. 
The hadronic system is reconstructed in 16 dif- 
ferent modes including one kaon and at most 3 
pions (and no more than 1 71" ). To discrimi- 
nate between all combinations in a given event, 
the combination with the lowest value of \ 2 — 
r m ES -m B \2 + / ae_ n2 j se l e cted. The rate for 

B — > r/X s events with p*, > 2 GeV/c is found to 
be (6.8±S;g ± 0.6±g;g) x 10~ 4 , where the first er- 
ror is statistical, the second is systematic, and the 
third corresponds to the uncertainty in the contri- 
bution from B — ► rj'D^ decays. This measure- 



ment confirms the large rate previously measured 
by CLEO. Such a large rate is generally believed 
to be due to the anomalous coupling of the 7/ to 
two gluons. 

We have also investigated B decays to three- 
body final states. Particularly promising is the 
decay B° — > 7r + 7r~7r°, which could provide a 
clean measurement of the CP-violating phase a 
via a Dalitz amplitude analysis. Due to the 
small expected signal yield, the current version 
of the analysis aims to determine the B de- 
cay rates to ^±(770)^^, p°(770)7r°, p± (U50)n* , 
/3°(1450)7r°, (charged scalar) ± 7r T , f°n°, and non- 
resonant 7T + 7r - 7r° final states. A significant signal 
is found only in the B° — > p ± (770)n T mode. 

Also promising for the measurement of a is the 
decay B° — > Oq (980)7r T , which is expected to be 
free of main tree contributions. The analysis re- 
lies on a neural network algorithm to select B° de- 
cays to a (980)^, with — > 777^ and 77 — * 77. 
A likelihood fit using the neural network output, 
77 mass, niES and AE extracts a signal yield of 
18. events with a corresponding statistical 
significance of 3.7 a. The 7771"^ invariant mass is 
not included in the fit because the parameters of 
the ao resonance are not well known. 

Finally, we also measured the branching frac- 
tion for B + decays into K*°ir + . Part of the in- 
terest in this process stems from the expectation 
that tree and penguin amplitudes have compara- 
ble magnitudes, which makes this decay an ideal 
candidate for the observation of direct CP viola- 
tion, see Eq. ||. The measured yields and branch- 
ing fractions for this and all other rare hadronic 
decays are summarized in Table ||. 

5. DIRECT CP VIOLATION 

We searched for CP violation in the decay of 
the flavor self-tagged modes listed in Table |[ The 
table shows the measured CP-asymmetries as de- 
fined in Eq. [|. For the B° — > K + ir~ decay mode, 
an on-resonance sample with an integrated lumi- 
nosity of 30.4 fb^ 1 is used In the case of 
the B° — > p ± 7r T mode, no tagging of the B° fla- 
vor is attempted and the asymmetry quoted in 
Table || corresponds to the asymmetry between 
the number of p + ir~ and p~ir + events. No evi- 




Table 3 

Summary of direct CP violation measurements. 
The first and second uncertainties are statistical 
and systematic, respectively. 

Mode_ A C p 

-0.044 ±0.076 ±0.012 
-0.07 ±0.08 ±0.04 
0.00 ±0.18 ±0.04 
-0.21 ±0.18 ±0.03 
-0.05 ±0.20 ±0.03 

-0.43±81o ± 0.06 
0.00 ±0.27 ±0.03 

-o.oi±8:i! ± 0.03 

-0.11 ±0.11 ±0.02 
-0.04 ±0.18 ±0.02 



dence for direct CP violation has been found but 
the precision is now better than 10% in some of 
the more abundant modes and, thus, these data 
become more useful in constraining models. 

6. SUMMARY AND OUTLOOK 

The large data sample accumulated at the 
SLAC PEP-II B Factory is already providing a 
substantial amount of interesting results. We are 
now sensitive to rare B decays with branching 
fractions as low as 4 x 10~ 6 . A search for di- 
rect CP violation in radiative penguin and rare 
hadronic decays finds no evidence for CP viola- 
tion in the decay process. It should be noted that 
some of the results presented here are preliminary. 

PEP-II is expected to deliver an integrated lu- 
minosity of approximately 100 fb _1 by the sum- 
mer of 2002, which will provide a five-fold in- 
crease in statistics compared to the results re- 
ported here. This will allow many more rare pro- 
cesses to be studied. 
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